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The  hardening  effects  of  focused  ion  beam  (FIB)  induced  damage  produced  during  the  fabrication  of 
micropillars  are  examined  by  introducing  a  surface  layer  of  nanosized  obstacles  into  a  dislocation  dynamics 
simulation.  The  influence  of  the  depth  and  strength  of  the  obstacles  as  a  function  of  pillar  diameter  is  assessed 
parametrically.  We  show  that  for  a  selected  set  of  sample  sizes  between  0.5  and  1.0  pm,  the  flow  strength  can 
increase  by  10-20  %,  for  an  obstacle  strength  of  750  MPa,  and  damage  depth  of  100  nm.  On  the  other  hand, 
for  sizes  larger  and  smaller  than  this  range,  the  effect  of  damage  is  negligible.  Results  show  that  the  obstacles 
formed  during  the  FIB  milling  may  be  expected  to  alter  the  microstructure  of  micropillars,  however,  they  have 
a  negligible  effect  on  the  observed  size-strength  scaling  laws. 
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Focused  ion  beam  (FIB)  methods  are  site  specific  sample 
preparation  techniques,  that  have  been  widely  used  for  a 
couple  of  decades  in  semiconductor,1'2  material  science, 3-9 
biology,1011  geoscience,12  and  many  other  disciplines.  In  FIB 
techniques,  the  shape  of  the  target  material  is  manipulated  by 
bombarding  the  surface  of  the  target  by  gallium  accelerated 
heavy  ions. 

One  possible  draw  back  of  the  FIB  process,  is  that  gal¬ 
lium  ions  can  critically  alter  the  defect  state  of  the  target 
material.  This  could  have  an  impact  on  samples  with  a  large 
surface-to-volume  ratio.  The  FIB  induced  damage  observed 
in  many  conventional  alloys,  can  be  present  in  the  form  of  an 
amorphous  layer,13,14  dislocation  networks,15,16  vacancy  clus¬ 
ters  in  the  form  of  staking  fault  tetrahedra,17  interstitial 
loops,  decoration  of  dislocations  with  Ga,  as  well  as  Ga 
precipitates.15,16  Moreover,  the  depth  and  concentration  of 
the  FIB  induced  damage  is  dependent  on  the  target  material, 
ion  acceleration  energy,  ion  dose,  and  the  ion  incident 
angle. 14,18-20 

In  a  recent  application  of  the  FIB  technique,  Uchic  et  al. 
introduced  a  new  experimental  technique  to  quantify  plastic¬ 
ity  size  dependency.21  Via  this  technique,  micropillars  were 
milled  from  a  bulk  crystal  using  a  FIB,  then  a  uniform  com¬ 
pression  load  was  imposed  on  the  top  of  the  micropillar  us¬ 
ing  a  nanoindentor  with  a  flat  ended  diamond  tip.  These  ex¬ 
periments  produced  an  unexpected  dependence  of  strength 
and  ductility  on  crystal  size.  This  sparked  a  wide  interest  in 
the  problem  of  plasticity  size  dependency  in  many  conven¬ 
tional  alloys.3-9  From  these  results,  it  was  suggested  that  the 
observed  size  effects  are  controlled  mainly  by  the  size  distri¬ 
bution  of  the  underlying  sources  and  dislocation 
microstructure.22 

Since  these  micropillars  have  a  large  surface-to-volume 
ratio,  the  effect  of  the  FIB  induced  damage  may  critically 
alter  the  mechanical  properties  of  these  pillars.  To  investigate 
this,  Kiener  et  al.,4  investigated  the  damaged  surface  due  to 
Ga  ion  bombardment  during  FIB  milling  of  Cu  micropillars. 
They  reported  a  substantial  depth  of  Ga  implantation  up  to 
50  pm.  They  also  performed  basic  theoretical  analysis  to 
estimate  the  influence  of  a  precipitation  hardened  damage 
layer.  It  was  finally  concluded  that  FIB  damage  limits  the 


applicability  of  microcompression  tests,  and  that  its  effect 
may  be  proportional  to  the  inverse  of  the  micropillar  size.  In 
addition.  Shim  et  al.9  have  shown  that  compressed  Mo-alloy 
pillars  in  the  as-grown  condition,  behave  like  dislocation  free 
crystals,  and  yield  at  near-theoretical  strength.  On  the  other 
hand,  micropillars  that  were  FIB  milled  before  compression, 
yielded  at  a  much  lower  strength.  These  studies  suggest  that 
the  FIB  altered  microstructure  could  play  a  significant  role  in 
the  micropillar  plasticity. 

On  the  other  hand,  Shan  et  al.6  performed  in  situ  TEM 
compression  experiments  on  Ni  single-crystal  micropillars. 
They  observed  that,  for  micropillars  having  sizes  below 
0.16  pm,  the  high  density  of  initial  defects  produced  by  FIB 
milling  were  completely  driven  out  of  the  crystal  by  the 
compression  loading.  The  near-theoretical  flow  strength 
computed  for  these  micropillars  was  observed  to  coincide 
with  the  elastic  behavior  of  source-free  mechanically  an¬ 
nealed  pillars.  This  suggests  that,  in  contradictory  to  the  the¬ 
oretical  estimates  of  Kiener  et  al.,4  the  precipitation  hardened 
damage  layer  has  a  minimal  effect  at  such  very  small  scales. 
These  experimental  studies,  as  well  as  studies  performed  on 
Au  micropillars,5  suggest  that  the  effects  of  the  FIB  milling 
is  negligible  if  the  original  crystal  has  some  initial  disloca¬ 
tion  structure  before  milling. 

To  study  the  size  effects  observed  in  micropillar  compres¬ 
sion  experiments,  a  number  of  three-dimensional  (3D)  dislo¬ 
cation  dynamics  (DD)  simulations  were  recently  performed. 
These  simulations  led  to  discovery  of  key  mechanisms  that 
appear  to  govern  the  experimentally  observed  size 
effects.23-25  Although,  the  results  of  such  simulations  are  in 
general  agreement  with  several  experimental  observations, 
all  such  DD  simulations  failed  to  address  the  effects  that  may 
arise  due  to  the  FIB  induced  damaged. 

To  address  this,  we  used  3D  DD  simulations  to  investigate 
the  effects  of  FIB  induced  damage  on  the  observed  micropil¬ 
lar  strength  size  dependence.  The  computational  method 
adopted  here  is  that  developed  in.26  In  this  formulation,  the 
boundary  element  method  is  coupled  with  the  3D  parametric 
dislocation  dynamics  formulation  to  incorporate  the  influ¬ 
ence  of  interfaces  on  dislocation  motion  and  to  describe  mi¬ 
croscopic  plastic  flow  in  finite  volumes.  The  simulations 
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mimic  the  compression  experiments  of  Dimiduk  et  al .,3  for 
compression  of  FCC  single-crystal  micropillars.  The  simu¬ 
lated  micropillars  are  oriented  in  the  [001]  direction,  have  a 
fixed  aspect  ratio  of  3:1,  and  the  pillar  diameters  are  in  the 
range  of  D  =  0.25-5.0  /Am.  The  initial  dislocation  density  in 
all  simulations  ranges  from  p=lX1012  to  15X1012  m-2, 
which  is  equivalent  to  those  reported  from  experiments.3-5  A 
compressive  stress  is  applied  to  the  top  surface  of  the  cylin¬ 
der,  while  the  bottom  surface  is  kept  fixed.  Mimicking  the 
loading  technique  of  the  experiments,  a  mixture  of  a  constant 
strain  rate  loading,  and  a  constant  stress  loading  is  applied. 
At  any  given  time  step,  if  the  calculated  plastic  strain  rate  is 
smaller  than  the  applied  strain  rate,  ep<e,  the  applied  stress 
is  increased  to  obey  &=E(e-ep).  On  the  other  hand,  if  the 
computed  plastic  strain  rate  is  greater  than  the  applied  strain 
rate,  sp>e,  the  applied  stress  is  held  constant,  d-=0.  Once 
this  condition  is  satisfied,  the  stress  is  kept  constant  until  the 
total  displacement  becomes  grater  than  the  applied  strain  rate 
multiplied  by  the  total  elapsed  time,  (e^H-s^,)  <  eAr.  Finally, 
the  applied  strain  rate  in  all  simulations  is  e  =  200  s-1 ,  and 
the  material  properties  chosen  here  are  those  of  Ni. 

Out  of  all  forms  of  FIB  induced  damage,  nanosized  Ga- 
based  precipitates  in  FCC  Cu  are  believed  to  have  the  largest 
effect  on  the  flow  strength.14  Thus,  the  current  analysis  will 
focus  on  the  presence  of  such  Ga  precipitates  in  the  damaged 
layer,  which  represents  the  worst  case  scenario.  These  pre¬ 
cipitates  are  treated  as  spherical  obstacles  having  a  radius 
rp=  1.4  nm.  Since  these  precipitates  have  a  much  smaller 
size  compared  to  the  volume  of  the  simulated  pillars,  the 
strain  field  of  these  precipitates  can  be  ignored.  In  addition, 
the  precipitates  are  randomly  distributed  in  the  outer  circum¬ 
ference  of  the  simulated  micropillar,  with  an  average  inter¬ 
precipitate  spacing  of  lp  =  10.2  nm,  and  within  a  layer  of 
depth  t  from  the  surface  of  the  pillar.  These  values  are  in 
agreement  with  experimental  observations  in  Cu  grains  that 
were  FIB  machined.15  The  layer  depth  is  also  varied,  such 
that  t  =  25,  50,  and  100  nm,  respectively.  These  values  are 
chosen  to  mimic  the  experimentally  predicted  depth,  which 
can  be  in  the  range  10^f^  100  nm,  depending  on  the  ion 
energy,  incident  angle  of  the  ions,  and  the  target 
material.14-18-20 

From  3D  dislocation  dynamics  simulations  performed  by 
Takahashi  and  Ghoniem,27  it  was  reported  that  the  critical 
resolved  shear  stress  (CRSS)  required  to  move  a  straight  dis¬ 
location  through  a  one-dimensional  precipitate  array  is 
mainly  dependent  on  the  size  of  the  precipitates,  interprecipi¬ 
tate  spacing,  and  the  ratio  of  the  precipitate-to-matrix  shear 
modulus.  From  their  results,  it  was  shown  that  the  CRSS  for 
spherical  precipitate,  having  diameters  in  the  range  2.5  <rp 
<5  nm,  can  reach  up  to  (0.3  —  0 .4)  X  /j,mb  I  lp  for  soft  pre¬ 
cipitates  having  /Ap//Am<S  1,  and  up  to  (0.6—0.85)  X  /imbl  lp 
for  strong  precipitates  having  /tp/ /tmS>  1.  On  the  other  hand, 
in  the  intermediate  range  where  0.1  < /A;,//Am<3,  the  CRSS 
is  in  the  range  (0  —  0.4)  X  /jimbllp.  Here,  fip  and  /a,„  are  the 
shear  moduli  of  the  precipitate  and  matrix,  respectively.  In 
addition,  the  precipitate  strength  was  shown  to  decrease  with 
decreasing  precipitate  size.  Since  the  size  of  the  precipitates 
in  our  simulations  are  much  smaller  than  those  used  in,27 
their  strength  would  be  expected  to  be  smaller.  Based  on 
these  results  and  due  to  the  lack  of  data  on  the  properties  of 


FIG.  1.  (Color  online)  Microstructure  of  a  D=0.15  /im  micro¬ 
pillar  at  0.17%  strain,  (a)  No  damage  layer,  (b)  Damage  layer  depth 
f=100  nm  and  precipitate  strength  of  ap=150  MPa  is  used.  The 
precipitates  are  not  plotted  for  clarity.  The  numbers  encircled  refer¬ 
ence  the  same  dislocations  in  both  simulations. 

the  nanosized  intermetallic  Ga  precipitates  forming  during 
FIB  milling,  we  chose  to  vary  the  strength  of  the  precipitates 
in  our  calculations  such  that  their  CRSS  is  in  the  range 
(0.1  — 0.4 )X  fimbl  lp.  Thus,  we  vary  the  CRSS  such  that  crp 
=  100,  250,  500,  and  750  MPa,  respectively.  When  a  disloca¬ 
tion  propagates  in  the  micropillar,  and  meets  one  of  these 
precipitates,  it  will  get  pinned.  The  dislocation  then  will  not 
cut  through  the  precipitate  until  the  local  resolved  shear 
stress  equals  the  CRSS  of  the  precipitate,  crp.  When  this  criti¬ 
cal  strength  is  reached,  the  dislocation  is  released  and  the 
precipitate  is  destroyed. 

Each  time  a  dislocation  shears  through  a  precipitate,  its 
strength  will  be  reduced.  Takahashi  and  Ghoniem27  reported 
that  the  strength  of  a  2.5-nm  precipitate  can  reduce  by  10% 
after  the  first  two  times  a  dislocation  shears  through  it.  Two 
more  passes  would  reduce  it  even  more  by  30%  of  its  origi¬ 
nal  strength.  Finally,  as  the  precipitate  is  completely  sheared 
through,  the  precipitate  strength  drops  by  40-45  %.  In  our 
current  simulations,  since  the  applied  load  is  not  allowed  to 
drop,  and  after  the  first  dislocation  shears  through  a  precipi¬ 
tate,  the  local  CRSS  at  this  precipitate  will  not  change  much 
for  subsequent  dislocations  arriving  at  the  location  of  this 
precipitate.  Also,  since  the  strength  of  the  precipitate  should 
be  reduced,  the  presence  of  an  already  sheared  precipitate  on 
subsequent  dislocation  motion  will  be  minimal.  Thus,  to 
minimize  the  amount  of  computations,  the  precipitate  can 
safely  be  assumed  to  be  destroyed  after  the  first  dislocation 
shears  through  it.  In  other  cases  were  the  applied  stress  is 
allowed  to  drop,  this  assumption  would  not  be  valid. 

Figure  1  shows  the  microstructure  after  0.17%  strain  for  a 
simulated  micropillar  having  size  £>  =  0.75  /Am.  The  initial 
dislocation  density  was  p=3.4X  1012  m-2.  Figure  1(a),  illus¬ 
trates  the  micropillar  simulated  without  a  FIB  damage  layer. 
In  this  case,  a  flow  strength  of  246  MPa  was  reached.  Figure 
1(b),  shows  the  same  simulated  micropillar  but  with  t 
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FIG.  2.  (Color  online)  Dislocation  density  versus  strain  for  a 
D  =  0.75  |itm  micropillar  having  two  precipitate  strengths:  (a) 
crp  =  500  MPa  and  (b)  crp=150  MPa.  Both  cases  produce  an  initial 
increase  in  dislocation  density,  however  this  density  persists  into  the 
steady-state  regime  for  the  larger  precipitate  strength. 

=  100  nm,  and  crp= 750  MPa.  A  higher  flow  strength  of  273 
MPa  was  reached.  In  this  case,  the  precipitates  acted  as  an 
additional  barrier  to  the  dislocation  motion  near  the  surface, 
and  thus  a  higher  strength  was  required  for  the  dislocation  to 
penetrate  through  this  damaged  layer.  Examples  of  observed 
microstructure  changes  from  Fig.  1  are:  (1)  sources  that  were 
activated  when  the  damage  layer  was  not  included,  are 
trapped  when  including  the  damage  layer  (dislocations  2  and 
3);  (2)  additional  sources  being  activated,  when  including  the 
damage  layer,  due  to  the  higher  flow  strength  levels  (dislo¬ 
cation  6);  (3)  some  dislocation  locks/dipoles  were  prevented 
due  to  the  changes  in  the  microstructure  (dislocation  4). 

In  Fig.  2,  the  dislocation  density  versus  strain  for  a 
£>  =  0.75  fim  micropillar,  with  an  initial  dislocation  density 
of  p= 6.4  X  1012  m-2,  is  shown.  Results  for  various  damage 
layer  depth,  and  crp= 500  and  750  MPa  are  shown,  respec¬ 
tively.  It  is  observed  that  for  crp=  750  MPa,  the  dislocation 
density  increases  by  a  factor  of  ~2.5  compared  to  the  simu¬ 
lation  without  the  damage  layer.  The  simulation  results  show 
a  sharp  increase,  then  decrease  to  a  steady  state,  in  the  dis¬ 
location  density  as  shown  in  Fig.  2.  As  the  largest  sources  in 
the  pillar  are  first  activated,  part  of  the  dislocations  will  get 
pinned  in  the  damage  layer,  while  the  remainder  of  these 
sources  are  free  to  expand  in  the  precipitate-free  region.  If 
the  local  strength  is  not  high  enough  for  the  dislocation  to 
penetrate  through  the  damage  layer,  a  dislocation  pileup 


FIG.  3.  (Color  online)  The  average  percentage  change  in  the 
computed  flow  strength,  for  different  micropillar  sizes  and  precipi¬ 
tate  layer  depth,  (a)  <7^=500  MPa;  (b)  ap=150  MPa.  The  error 
bars  are  the  standard  deviation  of  the  simulated  results  from  the 
average  and  the  solid  lines  represent  the  best  fit  for  each  damage 
depth. 

builds  up,  and  the  dislocation  density  sharply  increases.  After 
one  or  two  passes  from  those  activated  sources,  a  higher 
local  strength  is  reached  that  can  destroy  the  precipitates  and 
the  dislocation  pileups  escape  the  pillar.  Once  this  occurs,  the 
dislocation  density  rapidly  decreases  to  a  lower  steady-state 
value.  It  is  also  observed  that  some  activated  dislocations 
that  would  have  escaped  the  pillar  if  the  damage  layer  was 
removed,  will  instead  be  trapped  in  the  pillar.  The  number  of 
trapped  dislocations  increases  as  the  precipitate  strength  in¬ 
creases,  and  a  higher  steady-state  dislocation  density  is 
reached. 

Figure  3  shows  the  percentage  change  in  the  computed 
flow  strength  as  a  function  of  pillar  diameter,  for  crp=500 
and  750  MPa,  respectively.  Here,  the  percentage  change  is 
defined  as  [((jj^-ay)/oy]  X  100%,  where  cr^  is  the  com¬ 
puted  flow  strength  when  including  the  damage  layer,  and  oy 
is  the  flow  strength  for  the  same  pillar  and  initial  dislocation 
structure  but  with  no  damage  layer.  For  each  pillar  size  be¬ 
tween  0.25  and  1.0  /zm,  six  to  eight  different  initial  disloca¬ 
tion  distributions  with  different  initial  dislocation  densities 
were  simulated,  for  each  damage  depth  and  precipitate 
strength.  In  addition,  the  2.5  and  5.0  /jl m  pillars  were  simu¬ 
lated  with  three  and  two  different  initial  dislocation  distribu¬ 
tions  and  initial  dislocation  densities,  respectively. 

It  is  clear  from  Fig.  3  that  there  is  a  window  of  micropillar 
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FIG.  4.  (Color  online)  A  comparison  of  the  flow  strength  at 
0.5%  strain  versus  micropillar  size,  showing  results  from  simula¬ 
tions  with  no  damage  layer  (Ref.  25),  and  the  current  simulation 
results  with  a  layer  of  precipitates  having  strength  500  and  750 
MPa.  Results  for  all  three  layer  depths  are  included. 

sizes  for  which  the  flow  strength  is  dependent  on  the  damage 
layer.  For  micropillars  having  sizes  at  and  below  0.25  /rm, 
the  number  of  sources  in  the  pillar  are  limited,  and  the  ap¬ 
plied  strength  needed  to  activate  the  dislocation  sources  is 
already  high.2"’  Thus,  the  activated  dislocations  can  overcome 
the  precipitates  without  the  need  for  a  higher  strength.  In 
addition,  for  pillars  having  diameters  equal  to  or  greater  than 
2.5  yum,  the  number  of  sources  in  the  pillar  is  high,  and  the 
local  dislocation  stress  field  becomes  strong  enough,  due  to 
the  initial  dislocation  pileup,  to  overcome  the  damage  layer 
without  the  need  for  any  additional  increase  in  strength.  On 
the  other  hand,  for  simulated  pillars  having  sizes  between 
0.5- 1.0  /z m,  the  effect  of  the  damage  layer  is  observed  to 
increase  the  flow  strength  by  2.5-23  %,  depending  on  the 
depth  of  the  damage  layer,  and  the  strength  of  the  precipi¬ 
tates.  For  this  window  of  simulated  pillar  sizes,  the  compe¬ 
tition  between  the  numbers  of  sources  available,  local  dislo¬ 
cation  stress  field,  strength  required  to  activate  dislocation 
sources,  and  the  precipitate  strength,  all  comes  into  play. 
Simulations  with  a  precipitate  strength  of  100  and  250  MPa 
revealed  similar  observations,  but  the  increase  in  the  flow 
strength  was  limited  to  2.5-10  %  for  pillar  sizes  between 
0.5- 1.0  yttm.  Softening  was  also  observed  in  some  simu¬ 
lated  pillars  having  sizes  1.0  and  2.5  /zm.  This  is  because 
the  evolving  microstructure  (i.e.,  dislocation  speed  and 
shape)  changes  in  the  presence  of  the  damage  layer,  and 
dipoles/locks  that  would  have  formed  in  the  absence  of  the 
damage  layer,  are  instead  avoided.  Thus,  a  lower  strength 
would  be  required  to  activate  the  dislocation  system.  The 
maximum  softening  observed  was  2.5%. 

Figure  4  shows  a  comparison  of  the  computed  flow 
strength  versus  micropillar  size,  for  previous  DD  simulations 
that  do  not  include  a  damage  layer,25  and  the  current  results 
with  different  precipitate  strength  and  damage  layer  depth.  It 
is  observed  that  even  for  the  strongest  precipitate  strength 
and  damage  layer  depth,  the  strength  scaling  power-law  ex¬ 
ponent  is  at  best  weakly  dependent  on  the  presence  of  a 
damage  layer.  Including  the  damage  layer  may  alter  the  un¬ 
derlying  dislocation  microstructure,  however,  even  for  the 
strongest  precipitates  and  deepest  defect  layers,  this  produces 


less  than  a  25%  change  in  the  flow  strength.  This  is  an  order 
of  magnitude  smaller  than  the  effect  of  varying  the  source 
size  distribution,  which  can  lead  to  as  much  as  200-300  % 
change  in  the  flow  strength  for  each  pillar  size. 

Rao  et  al.23  showed  that  a  higher  initial  dislocation  den¬ 
sity  will  result  in  a  weaker  scaling  response  and  vice  versa. 
For  large  pillar  sizes,  simulations  having  a  higher  initial  dis¬ 
location  density  was  observed  to  result  in  a  higher  flow 
stress.  This  is  due  to  the  dominance  of  forest-hardening  pro¬ 
cesses  for  these  larger  pillars.  Such  trend  is  reversed  at 
smaller  sizes  were  higher  dislocation  densities  will  result  in  a 
weaker  pillar.  They  also  observed  that  the  scatter  is  consid¬ 
erably  less  for  high  dislocation  density  simulations  compared 
with  low-density  simulations.  Nevertheless,  the  observed 
scatter  is  still  an  order  of  magnitude  larger  than  the  effect  of 
the  damage  layer  reported  in  the  current  study.  Thus,  the 
dependency  of  the  flow  strength  on  the  initial  dislocation 
density  will  not  be  affected  by  the  presence  of  the  damage 
layer. 

It  should  be  mentioned  that  Kiener  et  al.4  performed  some 
simplified  estimates  based  on  the  precipitate  hardening 
knowledge  and  predicted  that  the  strength  of  the  pillars  can 
increase  significantly  over  50%  for  pillars  smaller  than 
1.0  /jl m.  This  is  in  fact  contradictory  to  all  experimental  evi¬ 
dence  that  show  that  for  such  small  pillars  the  FIB  induced 
damage  seems  to  have  a  relatively  small  effect.5,6  While  such 
simplified  estimates  may  provide  some  insight  on  the  role  of 
intermetallic  hardening,  they  defiantly  do  not  give  decisive 
conclusions  and  may  not  be  entirely  correct.  Such  simplified 
estimates  fail  to  take  into  account  the  complex  3D  disloca¬ 
tion  structures  that  are  present  in  the  pillar  and  the  complex 
stress  distribution  that  cannot  be  accounted  for  without  full 
3D  dislocation  dynamic  simulations.  The  results  of  the  cur¬ 
rent  study  are  observed  to  by  in  general  agreement  with  all 
experimental  observations  and  provide  some  insight  on  the 
mechanisms  involved  with  the  deformation  process. 

In  addition,  there  have  been  some  attempts  to  investigate 
the  effects  of  other  types  of  radiation  damage  on  the  multi¬ 
plication  and  motion  of  dislocations.  Rodney  and  Martin28 
performed  molecular-dynamics  simulations  to  investigate  the 
dislocation  pinning  by  interstitial  loops  in  nickel  crystals. 
They  concluded  that  the  dislocation  loops  represent  rather 
weak  dispersed  barriers  and  their  effect  on  hardening  was 
small.  Another  important  point  that  needs  thorough  investi¬ 
gation  is  the  effect  of  the  presence  of  an  amorphous  layer. 
This  layer  can  have  a  width  in  the  order  of  tens  of 
nanometer29  and  acts  as  a  barrier  to  dislocation  motion.  Sim¬ 
plified  two-dimensional  DD  calculations  were  performed  to 
study  the  effect  of  an  impenetrable  layer  on  the  motion  of 
dislocation  in  Cu  thin  films.30  It  was  shown  that  the  flow 
strength  increases  with  the  presence  of  an  impenetrable  layer 
by  about  20%  for  films  having  1.0  fim  thickness.  Although 
more  3D  DD  detailed  investigations  need  to  be  performed, 
such  results  would  suggest  that  the  effect  of  an  amorphous 
layer  would  be  the  dominant  effect  at  the  smaller  micropillar 
sizes. 

Finally,  more  detailed  experimental  investigations  are  still 
required  to  completely  characterize  the  exact  microstructure 
forming  due  to  the  FIB  milling.  This  would  aid  in  building 
and  performing  full  3D  DD  simulations  that  can  correctly 
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represent  the  experimental  microstructure,  which  may  result 
in  a  better  understanding  of  the  influence  and  the  relative 
importance  of  the  different  FIB  induced  defects  on  the  plas¬ 
ticity  of  micropillars. 

In  summary,  the  effect  on  plasticity  of  FIB  induced  dam¬ 
age  during  the  fabrication  process  of  micropillars  was  exam¬ 
ined.  A  layer  of  nanosized  precipitates  was  introduced  into 
dislocation  dynamic  simulations  to  represent  such  FIB  in¬ 
duced  damage.  The  influence  of  the  layer  depth  and  strength 
of  the  precipitate  was  then  analyzed.  It  should  be  noted  that, 
the  depth  of  the  damage  layer  is  a  function  of  the  incident 
angle  of  the  ion  beam,  the  ion  acceleration  energy,  and  the 
milled  material.  Thus,  the  current  results  should  be  viewed  as 
representing  a  qualitative  understanding  of  the  effects  of  the 
presence  of  such  a  layer  on  the  response  of  micropillars. 
We  show  that  for  a  window  of  simulated  micropillar  sizes 
between  0.5  and  1.0  /zm,  the  How  strength  can  increase  by 
10-23  %  for  cr;)=750  MPa  and  damage  depth  of  100  nm.  In 
addition,  the  dislocation  density  can  increase  by  a  factor  of 


~2.5.  For  pillars  sizes  that  are  outside  this  window 
(i.e.,  D  <0.5  /zm  and  D  >  1 .0  /rm),  a  negligible  effect  on 
the  flow  strength  is  observed.  Finally,  the  current  results 
show  that  the  obstacles  that  form  during  the  FIB  milling  do 
alter  the  microstructure  of  micropillars.  However,  even  for 
the  strongest  obstacles  and  deepest  defect  layers,  the  effect  is 
less  than  25%  of  the  flow  stress  expected  in  undamaged  mi¬ 
cropillars.  Since  the  observed  and  calculated  size  strengthen¬ 
ing  for  fee  metals  is  an  order  of  magnitude  larger  than  this,  it 
can  be  concluded  that  such  obstacles  that  could  form  in  these 
micropillars  have  a  relatively  small  effect  on  the  observed 
size  dependence. 
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